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Abstract: 

We perform a detailed analysis of the statistical properties of Poisson networks and show that 
the metric and topological properties of random cellular structures, can not be derived from 
simple models of random networks based on a poisson point distribution [1]. In particular we 
show that Lewis and Aboav-Wieare laws are not obeyed in these network. 



1 Introduction 

In the past two decades there has been a considerable interest in studying a class of non- 
equilibrium systems known collectively as Cellular Structures^!}. The geometrical and dynami- 
cal properties of these systems are best displayed in the familiar pattern formed by a soap froth 
confined between two transparent plates . Other examples include polycrystalline domains in 
metals, ceramics, magnetic domains and monolayers of fatty acids on surface of water. One 
can also mention other examples from material science, like cracking in glazes, fracture and 
dewetting of polymer films above their glass transition temperature [2]. Despite the diversity 
in systems in which cellular structures are formed, numerous experiments have shown that the 
long time statistical behaviour of these systems are charaterized by certain universal, system 
independent laws. This means that topological and geometrical constraints influence the prop- 
erties of these networks in a very essential way. For example, the fact that energy is associated 
with the length of the edges of the cells, immediately leads to the conclusion that in the two 
dimensional structures all the vertices are 3-valent. Combination of this result with the Euler 
character formula shows that the average number of sides per cell, (i) is equal to 6: 



Similar considerations in 3-dimensions, where energy is associated with the area of faces of cells, 
proves that all vertices are 4-valent and that : 
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where (/) and (i) are the average number of faces per cell and the average number of edges 
per face respectivly. Besides these properties which are a direct consequence of the Euler 
character formula, expriments have revealed a number of very general, properteis among which 
the most important are : i) von-Neumann [3], ii) Lewis [4] and hi) Aboav-Weare law [5,6]. Von- 
Neumann law refers to the rate of expansion of a single cell and is accounted for theoritically in 
a satisfactory way. It has also been generalized to curved surfaces where the curvature of the 
surfaces plays a role both in the dynamic of a single cell and in it's stability [7]. The Aboav 
Weare law which is statistical in nature and refers to the correlation of adjacent cells has been 
shown to hold in random 3-valent graphs, viewed as planer feynman graphs of 3 theory and 
solved by techniques of matrix models [8]. 

However there has been no explanation of the empirical Lewis law which states that the average 
area of cells has a linear relation with the number of cells, for large number of edges. 
However recently there have been some attempts [1,9] to derive both Lewis law and Aboav- 
weare law, from poisson networks, i.e. networks based on a poisson distribution of horizontal 
segments between a hxed set of parallel lines. It has been claimed [1] that in such networks 
both laws are obeyed. Also it has been claimed [1] that in one type of the 3-dimensional poison 
networks, 3-dimensional analogues of Lewis and Aboav-Weare law are obeyed. In this paper 
we study closely these poisson networks and show that the Lewis and Aboav-Weare law which 
is derived from these models contain highly non-linear terms, which is in sharp contrast with 
experimental data. 

The structure of this paper is as follows: in section 2, we introduce the poisson network in 
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general, such that we can use the result of this section for both 2-dimensional and 3-dimensional 
networks. In section 3 we show that Lewis and Aboav-weare laws are not obeyed in these net- 
works. In section 4 we show that even if one introduces an extra element of randomness into 
these models, namely, non-uniform density distribution, one can not obtain Lewis and Aboav- 
Weare law, although one can obtain better values for the moment of distributions of number of 
edges and faces in 2 and 3-dimensional cellular structure. Our coclusions are given in section 
5. 

2 Poisson networks with uniform density 

The two dimensional poisson networks are generated as follows (figure 1)[1]. One takes a family 
of parallel lines in the y-direction in the plane . The distance between lines does not affect the 
topological properties of the network, although it affectes the geometry . This distance is taken 
as uniform and equal to do- Suppose there are N columns , C±, C2, . . . , Cjv, in the network. 
The region between two successive columns is divided into cells. The division of each column 
C a ,is based on a poisson point distribution in the y-direction. An edge in the x-direction is 
taken through each Poisson point (P-point) . In this way a 3-valent network is generated in 
the plane which, although different from the realistic cellular structures, is simple enough to be 
analysed closely for the study of various statistical-topological propertis of the network. The 
3-dimentional tetravalent network are generated by a similar process (figure 2 )[1]. One takes 
an arbitrary 2-dimensional trivalent network (base network) on the xy plane (figure 2(a)) and 
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takes vertical planes (parallel to the z-axis) through each edge which divides the 3-dimensional 
space into prismatic columns. In each column one considers a Poisson distribution of uniform 
density in the z-direction and divides the columns into cells by planes perpendicular to the z- 
axis through each P-point (figure 2(b)). The height of a cell, L, is the distance between adjacent 
P-points in the associated distribution (figure 2(c)). The average cell height is unity. Consider 
one of the columns, say C a . The Poisson distribution of n points in a segment of length L is 
given by : 

p(n) = ^e~ L . (3) 

Cosider a particular set of columns C a , Cp, . . . , C 7 . The probability that in a given distance L, 
there are n a , ng, . . . , n 7 segments respectivly in columns C a ,Cp, . . . , C 7 , is equal to 

$(L, n a ,np, . . . , n 7 ) = P{n a )P(np) . . . P(n 7 ) (4) 

Now consider a reference column which we denote by C a . The probability that there is a cell 
with height between L and L + dL in this column is e~ L dL. We denote the columns which 
are neighbors of this reference column by C ai , C a2 , . . . , C ak . The number of these neighboring 
columns depend on the dimensionality and type of the network. The joint probability that 
there is a cell of height between L and L + dL in the column C a such that its neighbors 
C ai i C a2 , . . . , C ak respectively have n ai , n a2 , . . . , n ak points in the interval L is : 

)dL = P(n ai )P(n a2 ) . . . P(n ak )e- L dL (5) 

clearly this distribution function depends on the individual values of the variables n ai , n a2 , . . . , n c 
however in the sequel another probability distribution will be useful, that is ^>(L,I)dL were I 
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is defined as / = n ai + n a2 + . . . + n ak . The joint probability^(L, I)dL, is calculated as follows: 

* a (L, I)dL = *a(L, n ai , n a2 , . . . , n ak )dL = &^ e ~^ L dL. (6) 

where J2' means the sum over all n a _j's subject to the constraint their sum be equal to I. It's 
important to note that the distribution ty a (L,I) is obtained from a sum of distribution of the 
form ^ a (L,n ai ,n a2 , . . . ,n ak ) and is not equal to any of them. (compare with eq.(9a) of[l]). 

3 statistical properties of poisson networks 

3.1 Two dimensional networks 

In two dimensions, every column has two neighbors, hence in eq.(6), k = 2, and the number 
of edges of a cell (see fig.(l)) in a reference column is equal to : 

1 = 1+4: (7) 

Combining eq.(6) and (7) gives the probability distribution g(L,i) of finding i- cells of hight in 
the interval (L, L + dL): 

g(L,i)dL= ( [ e' ? ' L dL. (8) 
clearly this distribution function is normalized. 

POO 00 

/ £ 5 (L,i)dL = l. (9) 
From eq.(8) one can obtain the average height of i- cells: 

<L> i =j Q Lg(L,i)dL = (10) 



This equation, then can be used to calculate the average area of i- cells. For the poisson networks 
considered in [1], where the width of all cells are equal to do, we find: 

< A >,= do < L >,= ^^r 4 (11) 

Clearly this relation, due to it's non-linear term, is in sharp contrast with Lewis law. 
From eq.(8), one can obtain the total probability distribution of i-cells g(i): 

poo 1 O 

g(i) = J o g(L,i)dL = -{-)*-* (12) 
One can also define a generating function : 



oo 



G(q)=<e i «>='E i g(i)e i * (13) 

i=4 



A simple calculation shows that: 



e 4q 



G{q) = (14) 



or 



lnG(q) =4q- ln(3 - 2e q ) (15) 

From which one obtains by successive differentiation various connected moment of the ditribu- 
tion: 

<i>=6 (16) 
< i 2 > - < i >= 6 (17) 
< i 3 > -3 < i 2 >< i > +2 < % > 3 = 30 (18) 
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3.2. three-dimensional networks. 

Consider a 3-dimensional network and a column whose base is an i-cell. In eq.(6), we should 
now equate k with i and, relate I with the number of faces as follows: 

/ + 2 + i = F. (19) 

The probability distribution of F-cells (cells with F-faces) of height betweenL and L + dL, is 
now obtained as: 

m(L, F)dL = ( ^7-~2)! e ~ (m)L( * L (20) 
From which one can obtain the average height of F-cells whose base are i-cells. 

roc rr 7 1 j 

<L>,,= j o Lg ,(L,F )d L = 1 - TW (—f-^ (21) 
The average height of F-cells with any base is given by 

oo 

< L > F = J2 <L >F,i 9(i) (22) 
i=4 

where g(i) is the probability distribution of i-cells in the base. As an approximaion one can set 
i =< i >= 6 in (21) and obtain: 

<L> F =^( 6 -) F ~* (23) 

Compare with eq.(36) of [1]. The average volume of F-cells will than be: 

< v > F = S^l( 6 -)^ (24) 

where S is the average area of cells in the base. 
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Here again, it's seen that there is a pronounced non-linear term, which however is not as 
worse as the two dimensional case. From eq.(20) one can also find the total ditribution of F-cells 
with i-cell bases: 

^> = ^<^> F - 2 (25 > 

From this formula various monents can be calculated: 

< F >i= 2i + 2 (26) 

< F 2 >= 58 + hi 2 (27) 

and 

< F >= 2 < i > +2 = 14 (28) 

< F 2 >= 58 + 5 < i 2 >= 268 (29) 

In the remaining part of this section we review breifly the basic steps of the analysis of ref. [1] 
and show that if one uses the correct form of the average < L >j as given by eq.(10), then for 
a poisson network, Aboav-Weare law will not be obeyed either, even approximatly. 
Consider an i-cell a of length L in a two dimensional poisson network(figure 3). This cell is in 
column a and has two neighbors in this column, called a' and a". 

There are two adjacent neighbors b\ and &2j which respectively distribute n\ and n<i points 
inside the cell a. The number of sides of a is then equal to L 

i = ni + ri2 + 4 (30) 
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The total number of sides of the cells in column b± adjacent to a is 

Ji = 4(m + 1) + (mi + 1) + (mi + 1) + m'/ (31) 

where the meaning of the numbers m\,m! 1 and mi" are specified in fig. (3). Clearly irrii which 
is the total number of sides of the cells adjacent to a is: 

irrii = 12 + J x + h- (32) 

Where the definition of J 2 is similar to that of J\ . For the average value of mi , m\ and m'[ we 
use: 

< mi + m[ >= < Lc ~ La > =< L c -L a > (33) 

< m " >= < ^ > =< L Cb - L a > + < L a > (34) 

Where the segment L Cb in column b\, which is the smallest segment containig L a , is called the 
cvovering lenght of l a [1]. With the probability distribution found in [1] for L c it is shown that 
the average < L c — L a >= 2, it then follows that: 

< irrii >= 12 + 4(ni + 112)+ < mi + rn'i > + < m-2 + m! 2 > + < m'[ > + < m' 2 ' > +4 (35) 

The number 12 comes from the average of sides of a' and a". Combining eq.(33), (34) and (35) 
one finds that: 

< inn >= 16 + 4i + 2 < L >< (36) 
Combining eq.(36) with (10) one obtains 

< im t >= 16 + ±1 + ^^)*- 4 (37) 
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which shows that Aboav-Weare law is not valid for such poisson networks. Clearly Poisson net- 
works by construction do not have enough randomness to be good models of cellular structures, 
as far as geometrical propertis (i.e. Lewis law) are concerned. One may even try to introduce 
more randomness, (i.e. non- uniform point distribution in different columns) to see if the non- 
linear terms shown in eq.(37) can be smoothed out. In the next section we examine breifly 
this possibility. The result is that although one can use this extra degree of freedom to obtain 
better agreement with exprimental data for various moments of number of sides(i.e. < i 2 >, 
etc.), again the Lewis and Aboav-weare law will not be valid. 



4 Poisson networks with non-uniform density 

A little modification of the previous formulas is necessery to study the propertis of a network 
where the density of points (the average distance between points to be denoted by A, in the 
following) in different columns are different. Instead of (1) we will have: 

Pa (n) = (^) n (38) 

A a TV. 



Instead of (4) we will have: 

. . n,u\ = Pi (n.T^Pn(nn) . . . Pu(nu\ = 

a=l 



k L_ 1 _l_ 



ML,n 1 ,n2,...,n k ) = P 1 (n 1 )P2{n 2 )...P k (n k )=Y[(—) n -—e *« (39) 



defining ^(L,I)dL as before, we obtain: 

Vo(L,I)= ]T $(L,n 1 ,n 2 ,...,n k )—e--odL= { ° xJ e~ L ^ (40) 
ni+n2+ ... +nk=I A \ n 
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where £o = j^ + j^ + j^ + -- - + j^- in a two dimensional network every column say C n has 
two neighbors, C ra _i and C n+ \. The number of sides of i-cell(a cell with i-edges) in column n, 
is i = 4 + I, where I = n n _i + n n+ \ and n n _i + n n+ i are the number of vertices contributed 
by the cells in adjacent columns. Then for the distribution of i- cells of length between (L and 
L + dL) in column C n , we obtain: 

*>™- e (41) 

where £ n = t—^ h t- + t^- We will also obtain, the distribution of i- cells in column C n : 

r°° 1 1 

g n (i)= g n ( L ,i)dL = —(l -—)*-* (42) 

JO ^n^n ^rt?n 

Clearly this distribution function is normalized. We can now evaluate various moments of the 
distribiution g n (i). In all the following calculations, we use the following formulas for geometric 



sum: 

oo 



YV = — (43) 



fc=0 

oo 



£^<^p <44> 



A simple calculation shows that: 



oo x i 
< i >n=y2i9n(i) = 3 + \ n i n = 4 + A n (- + ) (46) 

j =4 A n+1 



where < i > n is the average number of i-cells in column n. For a uniform density, this number 
turns out to be equal to 6. Another quantity of interest is < i 2 > n = X^U ^9n if) , which turns 
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out to be : 

< i 2 > n = 2(A n £„) 2 + 5A„£ n + 9 (47) 
The variance turns out to be equal to: 

<i 2 > n - < i >l= {\nin)(Ktn ~ I) (48) 

In order to find the moments < % > and < i 2 > in the whole lattice and not in a particular 
column n, we replace the average over the cells in a network, by an average over the densities 
in an ensemble of networks. Thus from (46) and (47) we have : 

< i >= 4+ < A n (— — + ) > (49) 

<i 2 >=2< (X n Uf > +5 < > +9 (50) 
where the averages on the right hand side are performed with a suitable distribution 

P(\ 1 ,X 2 ,...,\ N ). 

For simplicity in the following we assume that this distribution, is a symmetric distribution, 
i.e: 

P(A 1 ,A 2 ,...,A iV ) = P(A 1 ,,A 2 ,,...,A J vO (51) 

where the primed indices are a permutation of the unprimed ones. Clearly this minor restriction 
allows our conclusions to be valid for a very large class of probability distributions. For a 
symmetric ditribution we can simplify (49) and obtain: 

< i >= 6 (52) 
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as expected. Note that for a uniform distribution, we obtain 

< i 2 >= 42 < i 2 > - < i > 2 = 6 (53) 

in agreement with [1]. From (41) we also obtain the average height of i-cells in the nth column: 

< L > nyt = / Lg n (L, i)dL = — (1 - —^y- 4 (54) 
Again one sees that the non-linear behaviour persists in this average. 



5 Three-Dimentional Networks 

Consider a 3-dimentional laminated poisson network, based on an arbitrary 2-dimentional net- 
work, and a perticular i-cell is the base. We label the column with this base by C^o- The 
number of faces in a cell in the column above this i-cell is 

F = i + 2 + J (55) 

where J = mi + 777,2 + . . . + m*, and is the number of additional lateral faces which results 
from the cells in the adjacent k— column (fig. 4). Let fi(J) be the fraction of the F cells in 
column Cifl. In formula (40), one should now replace k by i, and I by J to obtain: 

ML, f, i) = -j^jr 1 e " L?0 ( 56 ) 

where £0 = ^ + X7 + ij + -- - + a"- Here Ao is the density of P-point distribution in the column 
Cifl and A Q (a = 1, 2, . . . , i) are the density of P-point distribution in the adjacent columns. 
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^fo(L,F,i)dL is the probability of finding an F-cell in column C^o, whose height is between L 
and L + dL. Integrating over L, one obtains 

* (F,*) = -i-(l--i-) M (57) 

This is the probability of finding an F-cell in column C^o- One can now find various moments 
of this distribution^,: 

< F > i)0 = A £o + i + 1 (58) 

< F 2 > ij0 = (i + l) 2 + (2i + l)(A £o) + 2(A £o) 2 (59) 
Note that for uniform density one has :Ao£o = i + 1 and the above formula reduce to 

< F >i= 2i + 2 (60) 

<F 2 >i= (i + l)(5i + 4) (61) 

in agreement withQ. By the same averaging procedure as in the two dimensional case one can 
obtain 

< F >i= AoCo + i + 1 (62) 

< F 2 >,= (i + l) 2 + (2i + l)(A £o) + 2(A e ) 2 (63) 
The later averages are performed over the density distributions. 
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6 Conclusion 

Our main conclusion is that Lewis and Aboav-Weare laws do not apply to poisson networks. This 
is due to the rather rigid geometrical structure of these networks. We have shown that by- 
more randomizing these networks, (i.e. in the sence of random density distribution in different 
columns), one can not still obtain the above mentioned laws. In fact the only feature that 
these networks have in common with real cellular structures is that they are 3-valent, and the 
main property of real systems in which all the angles around any vertex are equal is oboviously 
absent in these networks.lt seems that the Lewis and to some extant Aboav-Weare law which 
are inherently geometrical in nature will not show up in those random models (i.e. Poisson) 
which neglect this very specific property of cellular structures. 
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